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Apparatus for Measuring Parameters of a Flowing Multiphase Mixture 



Cross-Reference To Related Patent Applications 

This application claims the benefit of U.S. Provisional Application No. 60/439,717 
(Cidra Docket No. CC-0552) filed January 13, 2003, which is incorporated herein by 
reference in its entirety. 

Technical Field 

This invention relates to an apparatus for measuring the parameters of a multiphase 
flow, and more particularly to an apparatus for measuring the speed of sound propagating in 
a multiphase mixture to determine parameters, such as mixture quality, particle size, 
vapor/mass ratio, liquid/vapor ratio, mass flow rate, enthalpy and volumetric flow rate of 
the flow in a pipe or unconfined space using ultrasonic transducers/sensors. 

Background Art 

In the exhaust of a LPT turbine, the quality of the steam has a significant impact on 
the efficiency of the turbine. If the steam is too dry, excess energy remains in the steam 
before it hits the condenser. If the steam is too wet, the excess liquid particles can damage 
the turbine and degrade the aerodynamics. Typically, the exhaust steam from an LPT 
turbine is nominally approx. 90% quality, or 10% wetness. It would be advantageous to 
provide a probe to enable the measurement of the steam quality of LPT turbines. 

The knowledge or determination of the different parameters of a process flow 
comprising a saturated vapor/liquid flow is used to provide feedback of the process to 
improve quality control of a process or detect problems or needed maintenance of the 
processing system. One such parameter of the vapor/liquid flow is vapor quality (e.g., 
steam quality) and "wetness" of the mixture. Vapor quality of a saturated vapor/liquid 
mixture is defined as the ratio of the mass of the vapor phase to the total mass of the 
mixture. Conversely, the "wetness" of a saturated vapor/liquid mixture is defined as the 
ratio of the mass of the liquid phase to the total mass of the mixture. 

Saturated mixtures exist at temperatures and pressures at which liquid and vapor 

phases coexist. The temperatures and pressures at which the liquid and vapor phases 

coexist lie under the "vapor bubble" (i.e., saturation lines) on a phase diagram. A 
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representative phase diagram for water is shown in Fig. 1. The collection of points known 
as the saturated liquid line and the collections of points known as the saturated vapor line 
define the vapor bubble. These two lines connect at, what is termed, the critical point. 
Saturated mixtures exist only under the vapor bubble. For pressures and temperatures 
5 outside of the vapor bubble, the fluid exists as a single phase and the properties of that fluid, 
such as density, enthalpy, internal energy, etc., are uniquely defined by the pressure and 
temperature. For common fluids, such as water, these properties are tabulated as functions 
of pressure and temperatures and are available through a variety of references including a 
website hosted by NIST (ref: http ://webbook.nist . go v/chemistrv/fluid/) . 

10 For fluids at pressures and temperatures that lie within the vapor bubble, the fluids 

represent mixtures of the liquid and vapor phase. Although the properties of both the vapor 
and liquid phases are well defined (and tabulated for known substances), the properties of 
the mixture are no longer uniquely defined as functions of pressure and temperature. In 
order to define the averaged properties of a saturated mixture, the ratio of the vapor and 

15 liquid components of the mixture must be defined. The quality of the mixture, in addition 
to the pressure and temperature, are defined and used to uniquely determine the properties 
of the mixture. 

Measuring the average properties of a mixture is important in many industrial 
application since it is the mass averaged properties of the working fluid that enter directly 

20 into monitoring the thermodynamic performance of many processes. For example, it is the 
difference in the flux of enthalpy of the steam mixture flowing into and exiting from a 
turbine that determines the maximum mechanical work that can be extracted from the 
working fluid, and thus is important to determine component efficiency. However, if the 
steam entering or exiting the turbine were saturated, pressure and temperature measurement 

25 would not be sufficient to determine the specific enthalpy, but rather, a measurement of the 
quality of the steam would be required to uniquely define the thermodynamic properties of 
the saturated steam mixture. Note that once the quality and pressure (or temperature) of a 
saturated mixture is defined, the thermodynamic properties of the mixture are defined 
through mixing laws provided the properties of the liquid and vapor sates are known, 

30 The present invention provides the means for measuring the speed of sound enables 

one to determine quality, which in turn enables one to calculate enthalpy, density, and other 
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properties of the mixture. In addition to measuring the specific enthalpy, a measurement of 
the total mass is also, in general, needed to determine the flux of enthalpy. 

There are many other situations where knowing the quality of a saturated mixture is 
beneficial. For example, in a steam power plant, the quality of the steam within the steam 
5 turbine affects blade life. Generally it is desired to operate so the quality is as high as 
possible throughout the turbine to minimize liquid water drops that will erode the metal 
blades. Knowing the quality at the turbine inlet and exhaust (or at the exhaust only if the 
inlet is super-heated) provides a means to monitor the quality throughout the turbine. Also, 
to monitor plant performance so that it can be operated at optimum conditions and to 

10 identify degradation effects, the steam turbine thermal performance must be known. This 
requires the fluid enthalpy at the inlet and exhaust of each turbine to be known. If the fluid 
at either or both locations is saturated, pressure and temperature measurements alone will 
not be enough to determine the enthalpy. However if an additional measurement of quality 
is made the enthalpy is then defined. In addition, there may be other applications in 

1 5 refrigeration cycles. 

The ability to measure the flow rate and composition of the saturated vapor / liquid 
mixtures within the conduits is an important aspect of any system or strategy design to 
optimize the performance of a system based on saturated vapor/liquid mixtures. The 
industry recognizes this, and has been developing a wide variety of technologies to perform 

20 this measurement. These include probe based devices, sampling devices, Venturis and 
ultrasonic devices 

Summary of the Invention 

Objects of the present invention include providing an apparatus for measuring the 
25 speed of sound propagating in a multiphase mixture to determine parameters of the flow in 
a confined (e.g. pipe, duct), unconfmed space or a pipe, for example, using ultrasonic 
transducers/sensors . 

According to the present invention, an apparatus for measuring at least one 
parameter of a mixture flowing through a pipe is provided. The apparatus includes a pair of 
30 ultrasonic transducers disposed axially along the pipe for measuring the transit time of an 
ultrasonic signal to propagate from one ultrasonic transducer to the other ultrasonic 
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transducer. A processor, responsive to said transit time signal, provides an output signal 
indicative of the at least one parameter of the mixture flowing through the pipe. 

The foregoing and other objects, features and advantages of the present invention 
will become more apparent in light of the following detailed description of exemplary 
5 embodiments thereof. 

Brief Description of the Drawings 

Fig. 1 is a representative phase diagram for water. 

Fig. 2 is a schematic illustration of a probe in accordance with the present invention. 
10 Fig. 3 is a perspective view of a probe in accordance with the present invention. 

Fig. 4a is a perspective view of a probe embodying the present invention mounted 
within a pipe having circular cross-section in accordance with the present invention. 

Fig. 4b is a perspective view of a probe embodying the present invention mounted 
within a duct having rectangular cross-section in accordance with the present invention. 
15 Fig. 5 is a cross-sectional view of a plurality of probes disposed within a pipe for 

characterizing the flow pattern of the flow passing through the pipe in accordance with the 
present invention. 

Fig. 6 is a side view of a steam turbine having a plurality of probes disposed at 
different stages of the turbine and a different depths within each stage in accordance with 
20 the present invention. 

Fig. 7 is a cross-sectional view of a probe in accordance with the present invention. 

Fig. 8 is a plot of sound speed as a function of frequency for vapor/liquid mixtures 
with fixed droplet size (50 mm) and varying vapor-to-liquid mass ratio in accordance with 
the present invention. 

25 Fig. 9 is a plot of sound speed as a function of frequency for vapor/liquid mixtures 

with varying particle size where the vapor-to-liquid mass ratio is equal to 1.8 in accordance 
with the present invention. 

Fig. 10 is a plot of sound speed as a function of frequency for vapor/liquid mixtures 
with varying particle size, in accordance with the present invention. 
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Fig. 1 1 is a flow diagram of an optimization procedure employed to determine 
vapor-to-liquid ratio and droplet size from analytical model and experimentally determined 
dispersive speed of sound data in accordance with the present invention. 

Fig. 12 is a plot of the speed of sound propagating through a saturated vapor/liquid 
5 mixture having varying temperature and pressures versus quality of the mixture, in 
accordance with the present invention. 

Fig. 13 is a plot of the volumetric vapor phase fraction for vapor/liquid mixtures 
having varying temperature and pressures versus quality of the mixture, in accordance with 
the present invention. 

10 Fig. 14 is a plot of the enthalpy/volume for vapor/liquid mixtures having varying 

temperature and pressures versus the speed of sound propagating through the mixture, in 
accordance with the present invention. 

Fig. 15 is a plot of the enthalpy/volume for vapor/liquid mixtures having varying 
temperature and pressures versus quality of the mixture, in accordance with the present 

15 invention. 

Fig. 16 is a schematic diagram of another embodiment of a probe embodying the 
present invention. 

Fig. 17 is a schematic illustration of another embodiment of a probe in accordance 
with the present invention. 
20 Fig. 18 is a graphical diagram of an ultrasonic signal scatters when contacting large 

particles in the flow in accordance with the present invention. 

Fig. 19 is a plot the acoustic wavelength versus the frequency in accordance with the 
present invention. 

25 Best Mode for Carrying Out the Invention 

Referring to Figs. 2 and 3, a probe, generally shown as 10, is provided to sense and 
determine specific characteristics or parameters of a multi-phase mixture 12 flowing 
through a pipe (conduit) or in an unconfined space. The multi-phase mixture may be a two- 
phase liquid/vapor mixture, a solid/vapor mixture or a solid/liquid mixture, or even a three- 
30 phase mixture provided the particles within the fluid are substantially small as will be 
described hereinafter. One example of a multiphase mixture that can be measured is a 
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saturated vapor/liquid mixture, such as steam. To simplify the description of the present 
invention, the probe 10 will be described as an apparatus for measuring the parameters of a 
steam mixture, however, one will appreciate that the probe may be used to measure specific 
characteristics of a particular multiphase mixture provide the particle disposed therein are 
5 small. 

As will be described in greater detail, the probe measures the speed of sound 
propagating through the multiphase mixture flow to determine any one of a plurality of 
parameters of the flow, such as the steam quality or "wetness", vapor/mass ratio, 
liquid/solid ratio, the volumetric flow rate, the mass flow rate, the size of the suspended 
10 particles, and the enthalpy of the flow. 

Fig. 2 illustrates a schematic drawing of the probe 10 that includes a sensing device 
16 comprising an ultrasonic transmitter 18 and an ultrasonic receiver 20 disposed axially 
along the pipe 14, spaced a predetermined distance. In response to a signal processor 22, 
the transmitter fires an ultrasonic signal 24 to the receiver which then provides a signal 26 
1 5 to the signal processor indicative of the transit time or time of flight of the ultrasonic signal. 
The signal processor 22, in response to the transit time measurement signal 26, provides a 
signal 28 to the speed of sound (SOS) processing unit 30. As will be explained in greater 
detail hereinafter, the transit time of the ultrasonic signal between the transmitter and the 
receiver is related to the speed of sound of the fluid 12. The SOS processing unit 30, in 
20 accordance with the relationships of the SOS to various parameters of the fluid flow, 
determines the desired parameters of the fluid flow, such as steam quality, particle size 
density and mass flow. Some of these relationships are illustrated in Figs. 8, 9, 10, and 12 - 
15. 

Typically, the majority of the mass of the liquid water in the LPT exhaust is 
25 contained in small water droplets, having a mean droplet diameter of -1/3 micron. These 
saturated steam mixtures can be classified as particulate flows in which small, low volume 
fraction of relatively incompressible, dense particles are contained in a continuous vapor 
phase. The particles affect the propagation of the sound through the mixture in a frequency 
dependent manner. Models of this propagation is described in US Patent Application Serial 
30 No. 10/349,716 (CC-0579), and US Patent Application Serial No. 10/ 376,427 (CC-0596), 
which are both incorporated herein by reference. In these models, the sound speed of the 
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mixture if a function of frequency and related to the volumetric phase fraction, viscousity 
and drag coefficient, and fluid density and particle density and sound speed, which will be 
described in greater detail hereinafter. 

Referring to Fig. 10, the frequency behavior can be categorized into three frequency 
5 regimes. In the low frequency regime, or quasi-steady regime, there is negligible slip 
between the particles and the gases during the propagation of the acoustic waves. In this 
regime, the liquid particles essential mass load the gaseous phase. In this quasi-steady 
regime the sound speed asymptotes to the quasi-steady value which is directly related to the 
quality of mixture. 

10 In the high frequency regimes, the particles essentially no longer participate in the 

propagation of the acoustic waves. The acoustic oscillations are too high frequency to 
accelerate and decelerate the particles, and the speed of propagation asymptotes to that of 
the gaseous phase in isolation. 

The quasi-steady and high frequency regimes are separated by a transitional regime 

15 in which the speed of propagation is highly dependent on frequency. Herein, the 

transitional frequency is defined as the frequency at which the speed of propagation is 
midway between the quasi steady propagation velocity and the high frequency propagation 
velocity. The transitional frequency is strongly affected by particle size. For a given steam 
mixture, the transitional frequency scales with the square of the inverse of particle size. 

20 The present invention utilizes the transit time method to determine the dispersive 

properties of steam mixtures to characterize steam quality and particle size. The invention 
further contemplates using the sing around method as shown in Fig. 17 to determine the 
dispersive properties. 

In addition to influencing the dispersive properties of the mixture, scattering of 

25 acoustic waves by the particles, which by impeding propagation of acoustics through the 
mixture, can adversely affect the ability of transit time devices to effectively measure the 
speed of sound of inhomogeneous mixtures. 

Acoustic methods, which rely on transmitting and receiving acoustic waves, can 
encounter difficulties with inhomogeneous flows due to scattering of the acoustics by the 

30 inhomogenieities. The scattering effect is shown schematically in Fig. 18. This typically 
occurs when the acoustic wavelength approaches the length scale (i.e., size) of the 
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inhomogenieities. If the wavelength of the acoustics is several orders of magnitude greater 
than the length scale of the acoustics, little scattering occurs as shown in Fig. 1. 

The wavelength of acoustics traveling in a mixture with a nominal sound speed of 
1200 fee/sec - representative of the speed of sound in slightly wet steam and the exhaust of 
5 an LP turbine is shown in Fig. 19. As shown, the wavelength of the sound in the range of 
10,000 Hz to 1,000,000 Hz is several orders of magnitude larger than the diameter of the 
steam particles for steam representative of conditions at the exhaust of an LPT. 

Based on the analysis herein, a transmitter and receiver capable of operating over a 
frequency range of 10,000 to 1,000,000 Hz would be capable of effectively measuring the 
i 0 dispersive properties of steam at the exhaust of LP for steam with particle size of 

approximately 1 micron and smaller. The same device would not be effective for steam 
with particle sizes 5 microns or larger. For mixtures with these particle sizes and larger, 
transmitters and receivers capable of operating at significantly lower frequencies would be 
required. 

15 For example, based on the dispersive properties given in Fig. 10, to effectively 

measure the dispersive properties of a steam mixture with particles (droplets) on the order 
of 50 microns, the acoustic transmitters and receivers would have to be capable of launched 
and deterring waves of frequency well below 100 Hz. 

The probe 10 may be used a number of different ways. For example as shown in 

20 Figs. 4a and 4b, the sensing device 16 of the probe may be mounted within a pipe 19 or duct 
23, respective to measure the fluid flow or mixture passing therethrough. The probe 10 is 
particularly useful for large diameter pipes or ducts having a large cross-sectional area, such 
as smokestacks, exhaust ducts or HVAC systems. The utility of the probe is especially 
evident for measuring the flow of a multiphase mixture 12 that is not confined within piping 

25 or ducting. For example, the probe may be mounted within a gas turbine to measure the 
steam "wetness" or other parameters of the steam exiting the exhaust duct of the steam or 
LPT turbine. 

The probe 10 is particularly useful in characterizing the pattern of the flow 12 within 
a confined space (e.g., a pipe 19, duct 21) or unconfined space. For example as shown in 
30 Fig. 6, one or more sensing devices 16 of respective probes may be disposed at different 
locations across the area of the flow 12. The data provided by each of the probes and the 
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known location of the probe enables one to characterize the flow pattern. For example, one 
probe may be disposed adjacent to the wall of a pipe 19 and another probe may be disposed 
at a central position within the pipe to characterize the velocity and vapor/mass (e.g. steam 
wetness) at the two locations. One will appreciate that any number of sensing devices 16 
5 may be disposed at number of locations within the flow. 

The invention also contemplates a probe having a plurality or array of 
interconnected sensing devices 16 at 25 that extends across the flow 12 in a single plane or 
disposed in different planes of the flow. In other words, the sensing devices 16 of the 
probes 10 may be disposed at different locations along the flow 12. For example as shown 

10 in Fig. 6, probes may be disposed at different stages of a gas steam turbine at different 

locations at each stage. The probe may be used to measure the steam wetness at each stage 
at different location or depth from the walls of the turbine within each stage to measure the 
efficiency of the turbine or aid with the design of the turbine. 

Referring to Figs. 2, 3 and 7, in one embodiment of the present invention, the 

15 sensing device 16 of the probe 10 includes an inner tube 14 disposed within a tubular, outer 
housing 50 to provide an input port 53 and output port 55 for the flow 12 to pass through. 
The inner tube is generally cylindrical in shape having an axial bore 57 with a circular 
cross-section. The invention, however, contemplates that the inner tube may be of any 
shape or cross-sectional shape, such as squares, oval rectangular or any other polygonal 

20 shape. The cross-sectional shape may even be different along the length of the inner tube. 
The outer diameter of the inner tube is approximately one inch, but the diameter may be of 
any length. 

A pair of end caps 52,54 is disposed at the respective ends of the inner tube 14 and 
outer housing 50 to maintain and support the tube coaxially within the housing. The housing 

25 protects the array of sensors 18-21 disposed along the tube from the flow 12, and also acts 
as an insulator or isolator to prevent external acoustic and/or unsteady pressure disturbances 
from affecting the sensors 18-21. The end caps have a central bore 56 with a inner diameter 
substantially the same as the inner diameter of the tube 14. The outer surface 58 of the end 
caps are tapered and the outer ends 60 are rounded to provide an aerodynamic profile to 

30 reduce the drag of the flow (e.g., steam) over the probe 10 to reduce the wind resistance and 
stresses thereon. The aerodynamic profile also reduces the disturbance of the flow of the 
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fluid or mixture. The aerodynamic characteristics are particularly important for high speed 
steam flow (e.g., 0.7 Mach), such as steam exiting a gas turbine exhaust. 

While the sensing device 16 has a pair of end caps 52,54 as shown best in Fig. 7, the 
invention completes a sensing device 16 having no discrete end caps and that the inner tube 
5 14 extends the length of sensing device of the probe 10. 

The embodiment of the sensing device 16 shows a single input and output port 
53,55, however, the invention contemplates that the sensing device may have a plurality of 
input ports and/or output ports that feed into and out of the central portion of the tube 14 
where the sensors array 18 - 21 are disposed. While contemplated, one will appreciate that 
10 additional drag may be place upon the sensing device 16 and additional disturbance to the 
flow 12. 

The sensing device 16 of the probe 10 further includes a fin-shaped support structure 
62 extending from the center of the housing 50 for mounting the sensing device to a wall or 
other support, as shown in Figs. 3 - 6. The support structure 62 includes a bore 64, disposed 

15 therethrough to communicate with the space 66 disposed between the tube 14 and the 

housing 50. The bore 64 provides a means to run the conductors of the ultrasonic sensors 
21 to the processing unit 24, as best shown in Fig. 2. The support structure 62 is oriented to 
reduce wind resistance to minimize disturbance of the steam flow. 

Steam driven turbines are a major source of electrical power world wide. At 

20 present, there are no real time, operationally effective methods to monitor the quality of the 
steam as it drives the generators. Ideally, the industry would like to extract the maximum 
amount of energy from the steam as it passes through the turbine system thus reducing it to 
water; however, as water droplets form from the steam, they induce erosion and wear in the 
turbine blades requiring expensive maintenance. Thus a balance between the energy 

25 extracted from and the amount of water entrained in the steam must be met. 

Standard temperature and pressure cannot uniquely determine the wetness, i.e. the 
amount of water in the steam, when both phases coexist. The present invention uses speed 
of sound measurements to determine wetness of saturated steam using dispersion 
calculations. In non-nuclear power generation plants, the steam is super heated, and thus 

30 comprises only one phase, for each process step except the Low Pressure (LP) Turbine 
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Exhaust. At this exhaust, the wet steam is traveling at speeds in the range 0.5 to 0.7 Mach 
(Ma). 

Referring to Fig. 2, a probe 10 embodying the present invention is provided that 
measures at least one parameter/characteristic of a multiphase mixture 12 such as a 
5 saturated vapor/liquid mixture 12 of liquid droplets suspended within a continuous 

vapor/gas, for example, flowing within a pipe 19, duct 21 or flowing unconfined (see Figs. 
3-6). The probe may be configured and programmed to measure the speed of sound 
propagating through the flow 12. Depending on the configuration or embodiment, the 
probe can measure at least one of the following parameters of the flow 12: the wetness or 

10 steam quality (volumetric phase fraction), the volumetric flow rate, the size of the liquid 

particles, the mass flow, the enthalpy and the velocity of the mixture. To determine any one 
of these parameters, the probe 10 measures the transit time of an ultrasonic signal through 
the fluid 12 between an ultrasonic transmitter 18 and receiver 20. 

In this case, the wavelength of the measured acoustic signal determines the sensor 

15 spacing. The desired wavelength of the measured acoustic signal is dependent upon the 
dispersion of particles in the mixture flow, which is dependent on the particle size, which 
will be described in greater detail hereinafter 

The probe 10 can be used in any application that carries liquid droplets suspended in 
a vapor/gas through a pipe, such as in paper/pulp, petroleum and power generation 

20 applications. For example, the present invention is well suited to measure the parameters 
(e.g. vapor/liquid ratio, particle size) for power generation systems. 

As one example, the present invention will be discussed in the context of a steam 
delivery system for power generation, but one will appreciate that the probe 10 can be 
applied to any number of other applications, as discussed hereinbefore. 

25 As described hereinbefore, the probe 10 of the present invention may be configured 

and programmed to measure and process the transit time of an ultrasonic signal between an 
ultrasonic transmitter 18 and receiver 20 propagating through the mixture to determine 
parameters of the mixture flow 12. One such probe 10 is shown in Fig. 3 that measures the 
speed of sound (SOS) of one-dimensional sound waves propagating through the 

30 vapor/liquid mixture to determine the composition the mixture, namely the "wetness" or 
steam quality of the mixture. The probe is also capable of determining the average size of 
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the droplets, velocity of the mixture, enthalpy, mass flow, steam quality or wetness, density, 
and the volumetric flow rate of the mixture. It is known that sound propagates through 
various mediums at various speeds in such fields as SONAR and RADAR fields. The speed 
of sound of a mixture within the inner tube 14 may be determined using a number of known 
5 techniques, such as those set forth in U.S. Patent Application Serial No. 09/344,094, entitled 
"Fluid Parameter Measurement in Pipes Using Acoustic Pressures'', filed June 25, 1999, 
and U.S. Patent Application Serial No. 10/007,749, entitled "Fluid Parameter Measurement 
in Pipes Using Acoustic Pressures", filed November 7, 2001, each of which are 
incorporated herein by reference. The present invention utilizes at least one probe 10 to 

10 determine various parameters of the saturated vapor/liquid mixture, wherein one of the 
parameters is the speed at which sound travels within in the flow, as will be more fully 
described herein below. 

In accordance with the present invention, the speed of sound propagating through 
the vapor/liquid mixture 12 is measured by passively listening to the flow with an array of 

15 unsteady pressure sensors to determine the speed at which one-dimensional compression 
waves propagate through a vapor/liquid mixture contained within the tube 14. 

For relatively well-mixed vapor / liquid mixtures in which the liquid phase exists as 
small droplets within a continuous gas phase, the flow can be termed mist flow. Assuming 
that the droplets of the vapor/liquid mixture are small enough and the acoustic frequencies 

20 and the frequencies of perturbations associated with the acoustics are low enough for the 
droplets of liquid to exhibit negligible slip (both steady and unsteady), the sound speed can 
be assumed to be substantially non-dispersive (that is constant with frequency) and the 
volumetric phase fraction of the mixture could be determined through the Wood equation: 

25 Pmix = £ + iPi 
i=i 

Pmix a lix .=■ Pfi] 
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For one-dimensional waves propagating within a vacuum backed tube 14 ( or a tube 
immersed in large volume of low impedance fluid such as air at atmospheric conditions) , 
5 the compliance introduced by the tube ( in this case a circular tube of modulus E, radius R 
and wall thickness t) reduces the measured sound speed from the infinite dimensional sound 
speed. The effect of the conduit is given by the following relationship: 



= — + cr where a = — 

Pmix C measured Pmix C mix Et 

10 

Utilizing the relations above, the speed at which sound travels within the 
representative vapor/liquid mixture is a function of vapor/liquid mass ratio. The effect of 
increasing liquid fraction, i.e. decreasing vapor/liquid ratio, is to decrease the sound speed. 
Physically, adding liquid droplets effectively mass loads the mixture, while not appreciably 
15 changing the compressibility of the air. Over the parameter range of interest, the relation 
between mixture sound speed and vapor/liquid ratio is well behaved and monatomic. 

While the calibration curves based on predictions from first principles are 
encouraging, using empirical data mapping from sound speed to vapor/liquid ratio may 
result in improved accuracy of the present invention to measure the vapor/liquid fractions of 
20 the mixture. 

The sound speed increases with increasing frequency and asymptotes toward a 
constant value. The sound speed asymptote at higher frequency is essentially the sound 
speed of air only with no influence of the suspended liquid droplets. Also, it is apparent 
that the sound speed of the vapor/liquid mixture has not reached the quasi-steady limit at the 
25 lowest frequency for which sound speed was measured. The sound speed is continuing to 
decrease at the lower frequency limit. An important discovery of the present invention is 
that the speed at which sound propagates through droplets suspended in a continuous vapor 
is said to be dispersive. As defined herein, the speed at which acoustic waves propagate 
through dispersive mixtures varies with frequency. 
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For accurately measuring sound speeds at ultra-low frequencies, the data suggests 
that utilizing a quasi-steady model to interpret the relationship between sound speed, 
measured at frequencies above those at which the quasi-steady model is applicable, and the 
liquid-to-vapor ratio would be problematic, and may, in fact, be impractical. Thus, the key 
5 to understanding and interpreting the composition of vapor/liquid mixtures through sound 
speed measurements lies in the dispersive characteristics of the vapor/liquid mixture. 

In accordance with the present invention the dispersive nature of the system utilizes 
a first principles model of the interaction between the vapor and liquid droplets. This model 
is viewed as being representative of a class of models that seek to account for dispersive 

10 effects. Other models could be used to account for dispersive effects without altering the 
intent of this disclosure (for example, see the paper titled "Viscous Attenuation of Acoustic 
Waves in Suspensions" by R.L. Gibson, Jr. and M.N. Toksoz), which is incorporated herein 
by reference. The model allows for slip between the local velocity of the continuous vapor 
phase and that of the droplets. The drag force on the droplets by the continuous vapor is 

15 modeled by a force proportional to the difference between the local vapor velocity and that 
of the liquid droplets and is balanced by inertial force: 



where K = proportionality constant, U f = fluid velocity, U p = liquid droplet velocity, p p = 
20 liquid droplet density and v p = particle volume. 

The effect of the force on the continuous vapor phase by the liquid droplets is modeled as a 
force term in the axial momentum equation. The axial momentum equation for a control 
volume of area A and length Ax is given by: 



where P = pressure at locations x and Ax, (j) p = volume fraction of the liquid droplets, pr = 
vapor density. 



F dng =K(U f -U p ) = p p v p 



dt 
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The droplet drag force is given by: 



F*. = W f ~U P ) = C d A p X -p f (u f -U p } 



5 where C d = drag coefficient, A p = frontal area of liquid droplet and pf = vapor density. 

Using Stokes law for drag on a sphere at low Reynold's number gives the drag coefficient 

as: 



Re p f {U f -U p )D p 



15 



where D p = droplet diameter and p = vapor viscosity. 



Solving for K in this model yields: 



K = lnfjD p 



Using the above relations and 1 -dimensional acoustic modeling techniques, the following 
relation can be derived for the dispersive behavior of an idealized vapor/liquid mixture. 



20 



1 + - 



<PpPp 



Pf 



2..2 \ 



1 + 6) 



2 Pp V P 

K 2 



In the above relation, the fluid SOS, density (p) and viscosity (<f>) are those of the pure phase 
25 fluid, v p is the volume of individual droplets and ty p is the volumetric phase fraction of the 
droplets in the mixture. 
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Two parameters of primary interest in steam measurements are droplet size and 
liquid-to vapor mass ratio (i.e., steam quality or steam wetness). To this end, it is of interest 
to examine the dispersive characteristics of the mixture as a function of these two variables. 
Figs. 8 and 9 show the dispersive behavior for vapor/liquid mixtures with parameters typical 
5 of those used in steam flow systems. 

In particular Fig. 8 shows the predicted behavior for nominally 50 ^m size liquid 
droplets in vapor for a range of liquid-to- vapor ratios. As shown, the effect of liquid-to- 
vapor ratio is well defined in the low frequency limit. However, the effect of the liquid-to- 
vapor ratio becomes indistinguishable at higher frequencies, approaching the sound speed of 

10 the pure air at high frequencies (above -100 Hz). 

Similarly, Fig. 9 shows the predicted behavior for a vapor/liquid mixture with a 
liquid-to-vapor ratio of 1.8 with varying liquid droplet size. This figure illustrates that 
liquid droplet size has no influence on either the low frequency limit (quasi-steady ) sound 
speed, or on the high frequency limit of the sound speed. However, droplet size does have a 

1 5 pronounced effect in the transition region. 

Figs. 8 and 9 illustrate an important aspect of the present invention. Namely, that 
the dispersive properties of mixtures of droplets suspended in a continuous vapor can be 
broadly classified into three frequency regimes: low frequency range, high frequency range 
and a transitional frequency range. Although the effect of droplet size and liquid-to-vapor 

20 ratio are inter-related, the predominant effect of liquid-to-vapor ratio is to determine the low 
frequency limit of the sound speed to be measured and the predominate effect of droplet 
size is to determine the frequency range of the transitional regions. As droplet size 
increases, the frequency at which the dispersive properties appear decreases. For typical 
steam applications, this transitional region begins at fairly low frequencies, ~ 2Hz for 50 |im 

25 size particles. 

In the low frequency regime, the liquid droplets exhibit negligible slip with the 
vapor. The frequency range for which the no-slip, quasi-steady approximation is valid is a 
function of a variety of parameters including droplet size, continuous phase viscosity, 
droplet shape and droplet density. 
30 The quasi-steady sound speed is given by the low frequency limit of the above 

relation, where VLR is vapor/liquid ratio: 
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Note that droplet size does not affect the low frequency limit of the sound speed. 

Similar to Fig. 9, Fig. 10 shows the predicted behavior for a saturated vapor/liquid 
mixture with a liquid-to-vapor ratio with varying droplet size. Specifically, the particle 
sizes of three different mixtures include 50 urn, 20 urn and lum. The transitional frequency 
range of the mixture having 50 urn droplets is approximately 3- 13 Hz, a central frequency 
(fi/ 2 ) of approximately 8 Hz. The transitional frequency range of the mixture having 20um 
droplets is approximately 1 1 - 1 10 Hz, a central frequency (f l/2 ) of approximately 60 Hz. 
The transitional frequency range of the mixture having 1 urn is approximately 8 - 80 KHz, a 
central frequency {f m ) of approximately 40 degrees. As shown, the droplet size greatly 
influences the dispersion characteristics of the saturated vapor/liquid mixture. The 
transistion from the quasi-steady state to the high frequency regime scales inversely with 
the square of the droplet diameter. As discussed hereinbefore, the dispersion characteristics 
set the frequency requirements for measuring the speed of sound propagating through the 
mixture to measure parameters of the mixture. 

The frequency of the speed of sound that is detected for a particular mixture sets the 
wavelength of interest. The wavelength is the inverse of the frequency, and therefore, the 
higher the frequency, the shorter the wavelength and vice versa. As described hereinbefore, 
the aperture should be at least a significant fraction of the length of the wavelength of the 
speed of sound of interest. A vapor/liquid mixture having droplets of approximately 3 urn 
has a central frequency (f 1/2 ) of approximately 3 KHz. 

In the high frequency limit, the dispersion relation predicts the sound speed with 
asymptote towards the sound speed of the pure vapor. 

=> «>) = <* fluid 

Interestingly, the high frequency limit is independent of both droplet size and liquid-to- 
vapor ratio. 
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Given the difficulties measuring sufficiently low frequencies to apply the quasi- 
steady model and recognizing that the high frequency sound speed contains no direct 
information on either droplet size or liquid-to-vapor ratio, it becomes apparent that the 
dispersive characteristics of the vapor/liquid mixture should be utilized to determine droplet 
5 size and liquid-to- vapor ratio based on speed of sound measurements. 

As described hereinbefore, the probe 10 of the present invention includes the ability 
to accurately determine the average droplet size and the liquid/vapor ratio of the liquid in 
the vapor/liquid mixture. Provided there is no appreciable slip between the vapor and the 
liquid droplet, the propagation of one-dimensional sound wave through multiphase mixtures 

10 is influenced by the effective mass and the effective compressibility of the mixture. For an 
air transport system, the degree to which the no-slip assumption applies is a strong function 
of droplet size and frequency. In the limit of small droplets and low frequency, the no-slip 
assumption is valid. As the size of the droplet increases and the frequency of the sound 
waves increase, the non-slip assumption becomes increasing less valid. For a given average 

15 liquid droplet size, the increase in slip with frequency causes dispersion, or, in other words, 
the sound speed of the mixture to change with frequency. With appropriate calibration the 
dispersive characteristic of a mixture will provide a measurement of the average droplet 
size, as well as, the vapor to liquid ratio of the mixture. 



20 experimentally determined sound speed as function of frequency, the present invention 
includes an optimization procedure to simultaneously determine droplet size and VLR in 
liquid/vapor mixtures: 



Referring to Fig. 1 1 there is shown an optimization procedure in accordance with the 
present invention in which the free parameters of an analytical model are optimized to 
minimize an error function. For illustration purposes, the error function utilized is the sum 



Using the model described above, which yields the equation shown below, and 




25 
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of the differences of the sound speeds between an analytical model and the experimentally 
determined sound speed as a function of frequency: 

err= £ W)^, -a(f) mt&smed ) 2 

f=fhw 

5 

Thus, the sound speed of a two-phase mixture varies with the ratio vapor and liquid 
phases present in the mixture. Through these relations, and using tabulated values for the 
sound speed and densities of the liquid and vapor phases of a process mixture, one can 
construct an explicit relationship between mixture sound speed and mixture quality. It 

10 should be noted that the Wood equation is an engineering approximation, the accuracy of 
which is dependent on the validity of a variety of assumptions. Experimental data may be 
required to define between quality and sound speed within required, but to be defined, 
accuracy limits. Various curves are produced in Fig. 12 showing the relationship of sound 
speed versus steam quality for well-mixed saturated steam mixtures over of range of 

1 5 temperatures and pressures. 

As is known in the art, the relationship between quality of a vapor / liquid mixture, a 
mass ratio, and the volumetric phase fraction of the vapor phase is dependent on the 
properties of the vapor and liquid phases. For steam the relationship is shown in Figs. 12 
and 13. According to an empirical flow model, the assumption of well mixed, mist-like 

20 flows are typically applicable for process mixtures having vapor volumetric phase fractions 
greater than 0.83 and with mixture velocities exceeding 3.5 * sqrt( D*g), where D is the 
tube 14 diameter and g is the acceleration due to gravity. For example, an 18 inch diameter 
steam tube translates to mixture velocities greater than -8 m/s (-26 ft/sec). 

As developed above, determining the enthalpy flux of a steam mixture is an 

25 important measurement. In accordance with the present invention when the total volumetric 
flow of the mixture is known, the enthalpy per unit volume of the mixture is needed to 
determine the total flow rate. Fig. 14 shows the relationship between enthalpy per unit 
volume and mixture sound speed for steam at a variety of conditions. The present invention 
further utilizes the relationship between enthalpy per unit volume and mixture sound speed 
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for steam at a variety of conditions and the relationship between enthalpy per unit volume 
and steam quality as shown in Fig. 15 to determine the quality of steam of a flow. 

In addition to measuring the liquid to vapor ratio of the mixture 12 and droplet size 
of the liquid suspended in the mixture using the measured speed of sound, the probe 10 
5 further includes the ability to measure of volumetric flow rate of the mixture by comparing 
the difference of the speed of an ultrasonic signal propagating with the mean flow and an 
ultrasonic signal propagating against the mean flow. 

The interaction results in sound waves propagating with the mean flow traveling at 
the speed of sound (if the vapor/liquid mixture were not flowing) plus the convection 
10 velocity and, conversely, sound waves traveling against the mean flow propagating at the 
speed of sound minus the convection velocity as illustrated in Fig. 2. That is, 

clr = a mix + u 
&l = amix - u 

15 

where a R = velocity of a right traveling ultrasonic signal relative to a stationary observer 
(i.e. the tube 14), a L = velocity of a left traveling ultrasonic signal apparent to a stationary 
observer, a m j x = speed of sound traveling through the mixture (if the mixture was not 
flowing) and u = the mean flow velocity (assumed to be flowing from left to right in this 
20 instance). Combining these two equations yields an equation for the mean velocity, 




Therefore, by measuring the propagation velocity of ultrasonic signals in both directions 
25 relative to the stationary tube 14 as described hereinbefore, the mean flow velocity can be 

calculated by multiplying the mean flow velocity by the cross-sectional area of the tube 14 
The practicality of using this method to determine the mean flow is predicated on 

the ability to resolve the sound speed in both directions with sufficient accuracy to 

determine the volumetric flow. For typical vapor liquid measurements, flow velocities are 
30 typically at ~ 10 ft/sec and sound speeds of -4000 ft/sec. Thus axial mach numbers are on 
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the order of 10/4000 of 0.0025. For a +/- 10% accuracy in flow rate (+/- 1 ft/sec), the sound 
speed of the upstream and downstream propagating waves would need to be resolved to +/- 
0.5/4000 or 1 part in 8,000. 

Fig. 17 illustrates another method of determining the transit time measurement know 
5 as the sing around method. As shown a series of ultrasonic transmitters and receivers 
propagate a series of ultra sonic signals wherein the receivers measure the transit time. 

While a couple of methods have been described to measure transit time of an 
ultrasonic signal through a fluid, the present invention contemplates that any method of 
measuring the transit time of an ultrasonic signal may be used to determine the speed of 
1 0 sound of the fluid for determining parameters of the fluid. 

However, for saturated vapor/liquid mixture flows, axial flow velocities are 
nominally around 70 ft/sec with no flow sound speeds of -700 ft/sec. This results in mach 
numbers of -0.1, approximately 2 orders of magnitude greater than typical vapor flows. For 
saturated vapor/liquid flows, to resolve the flow rate to 10% accuracy (or +/- 7 ft/sec), one 
15 would have to resolve the sound speed to +/- 3.5 ft/sec, or 3.5/700 or 1 part in 200. 

While the apparatus 10 of Figs. 2 - 7 has been described as a probe, the present 
invention contemplates that the apparatus may be use to measure fluid flow within a pipe, 
such as an industrial process pipe. 

It should be understood that any of the features, characteristics, alternatives or 
20 modifications described regarding a particular embodiment herein may also be applied, 
used, or incorporated with any other embodiment described herein. 

Although the invention has been described and illustrated with respect to exemplary 
embodiments thereof, the foregoing and various other additions and omissions may be made 
therein and thereto without departing from the spirit and scope of the present invention. 

25 
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